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STRACT 


A  nonlinear  damage  relation,  containing  the  axial  strain  history 
and  a  time  integral  over  the  stress  history,  is  proposed  for  the  case 
of  one-dimensional  time  dependent  tensile  stress.  Nonlinear  steady  and 
transient  creep  terms  are  included  in  the  axial  strain  relation,  and 
elastic  and  creep  Poisson's  ratios  are  introduced  into  the  lateral 
strain  relation.  Using  these  relations,  complete  damage  solutions  are 
obtained  for  the  constant  stress  rate,  step  stress,  relaxation  and 
constant  load  tests.  Observations  are  made  concerning  the  associated 
rupture  times. 
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1.  INTRODUCTION 


The  problem  of  creep  damage  and  rupture  is  currently  of  great 
practical  importance,  especially  in  the  field  of  nuclear  safety. 

This  subject  is  very  complex  and  much  remains  to  be  accomplished  in 
the  development  of  a  complete  and  useful  creep  mechanics  model  for 
the  estimation  of  rupture  time.  The  publication  based  on  the  Nov. 

1978  Ispra  course  on  creep  [1]  contains  detailed  discussion  on  every 
aspect  of  this  subject  as  well  as  an  exhaustive  bibliography;  also 
see  Hult  [2]  and  the  books  by  llult  [3]  and  Odqvist  [4].  In  a  recent 
contribution  to  the  subject  of  creep  rupture,  Belloni,  Bernasconi  and 
Piatti  [5]  have  employed  material  density  variation  as  a  measure  of 
damage  in  steel  under  constant  one-dimensional  tensile  load.  By 
extensive  data  analysis  they  found  that  the  damage  depends  on  the 
level  of  the  initial  stress,  the  time  explicitly  and  also  on  the  axial 
creep  strain  as  a  function  of  time.  By  employing  a  damage  relation 
expressed  in  terms  of  power  functions  in  these  three  variables  and  an 
axial  creep  strain  relation  consisting  simply  of  a  stress  power  law 
for  steady  creep,  they  were  able  to  estimate  the  creep  rupture  time 
and  thereby  assign  a  physical  interpretation  (i.e.,  in  terms  of  density 
variation)  to  the  Kachanov  rupture  parameter  (see  [6]).  Further  results 
on  this  approach  to  creep  damage  were  presented  at  two  recent  congresses 
17,8]. 

In  this  paper  we  discuss  in  detail  an  extension  c{f  the  damage 
relation  presented  in  [5,7,8]  to  the  case  of  time  dependent  one¬ 
dimensional  tensile  stress.  We  first  modify  this  relation  by  stipulat¬ 
ing  that  it  be  employed  for  constant  stress  rather  than  constant  load, 
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and  then  postulate  that  for  variable  stress  the  damage  rate  be  a 
function  of  the  axial  creep  strain  rate,  axial  creep  strain,  damage 
and  stress.  Also,  we  employ  a  fairly  general  axial  strain  relation 
containing  linear  elastic,  nonlinear  transient  creep  and  nonlinear 
steady  creep  terms.  The  temperature  is  assumed  to  be  constant  in 
both  the  damage  and  strain  constitutive  relations.  Using  these 
relations,  damage  solution  are  then  obtained  for  several  prescribed 
one-dimensional  stress  and  strain  histories  as  well  as  for  the  case  of 
constant  load  accompanied  by  lateral  contraction.  Such  solutions  are 
useful  in  the  design  and  analysis  of  tensile  test  programs,  and  also  for 
the  estimation  of  rupture  times. 

Section  2  contains  the  one-dimensional  model  for  damage  and 
strain  (both  axial  and  lateral) ,  first  for  constant  stress  and  then 
as  extended  to  variable  stress.  In  Section  3  we  obtain  complete  damage 
solutions  for  the  constant  stress  rate  test,  the  step  stress  test  and 
the  relaxation  test  with  transient  creep  excluded  in  this  last  case. 
Solutions  are  then  obtained  in  Section  A  for  the  constant  load  test 
with  lateral  contraction,  first  for  large  deformation  with  transient 
creep  excluded  (or  included  approximately  in  the  elastic  term)  and  a 
second  for  small  deformation  including  transient  creep  as  a  separate 
term.  A  short  discussion  containing  some  suggestions  for  possible 
testing  programs  then  concludes  the  paper. 
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2.  CREEP  MODEL 

2.1  Damage  Relation 

For  constant  one-dimensional  tensile  stress  at  constant  tempera¬ 
ture  we  employ  the  damage  relation 


D(t)  =  C  ec(t)a  oj  t'S 


(1) 


Here,  D(t)  and  ^(t)  are  t^ie  damage  and  axial  creep  strain  respec¬ 
tively  at  time  t  ,  is  the  constant  stress,  and  C,a,y,6  are 

positive  material  damage  constants  at  a  particular  temperature.  In 
[5]  was  the  initial  stress  for  constant  load  and  the  damage  D(t) 

was  given  as  -Ap/p^  where  Ap  is  the  density  change  p ( t )  -  pQ 
(negative)  due  to  void  formation  and  growth;  here  D(t)  is  not  necessar¬ 
ily  restricted  to  this  particular  measure  of  damage.  Equation  (1) 
may  also  be  written  as 


d_ 

dt 


D(t) 


C  ec(t) 


a 


(2) 


Note  that  cc(t)  does  not  include  the  elastic  strain,  which  is  con¬ 
sistent  with  the  view  that  damage  (such  as  due  to  void  formation  and 
growth)  is  not  associated  with  recoverable  deformation. 

In  generalizing  eq.  (1)  to  the  case  of  variable  stress  o(t),  we 
postulate  that  the  damage  rate  be  a  function  of  the  axial  creep  strain 
rate,  axial  creep  strain,  damage  and  stress,  i.e. 


D(t)  =  f[cc(t),  cc(t),  D(t),  o(t)l 


MWWif  U  ./• 


(3) 
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where  the  dot  above  a  variable  indicates  d/dt.  An  expression  consis¬ 
tent  with  this  function  is  obtained  directly  if  we  simply  replace  oq 
in  eq.  (2)  by  o(t)  ,  which  upon  integration  yields  the  desired  result 


D(t)  =  C  e c ( t ) 01 


a(t')Y/,Sdt' 


(4) 


For  a(t)  =  oq  eq.  (4)  clearly  simplifies  to  eq.  (1),  and  furthermore 
it  possesses  the  physically  necessary  characteristic  that  D(t)  be 
continuous  even  if  the  stress  history  is  discontinuous. 


2 . 2  Strain  Relations 

For  the  axial  creep  strain  at  constant  tensile  stress  and  temper¬ 
ature  appearing  in  Eq.  (1),  we  employ  in  this  paper  the  sum  of  axial 
creep  strain  components 


£  (t)  =  c  ( t )  +  £  (t) 
c  s  t 


where 


rtU) 


) 


(5) 


(6) 


Here,  Eg(t)  and  t  (t)  represent  the  axial  strains  due  to  steady 
creep  and  transient  creep  respectively,  and  X,n,u,q,T  are  positive 
material  creep  constants.  Expressions  of  this  type  have  been  widely 
used  in  the  study  of  creep  in  various  metals  and  polymers,  including 
the  practically  important  case  of  stainless  steel  (see  Garofalo  [9]) 


V. 
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The  power  form  in  time  is  also  widely  used,  but  for  the  sake  of 
brevity  we  restrict  our  attention  here  to  strain  in  the  form  of  eqs. 
(5-6). 

In  generalizing  eqs.  (6)  to  variable  stress  o(t)  we  follow 
the  strain-hardening  approach  of  Cozzarelli  and  Shaw  [10],  whereby 
each  component  of  creep  strain  rate  is  expressed  as  a  function  of  the 
same  component  of  creep  strain  and  the  stress,  i.e. 


f- s ( t )  =  g[eg(t),  o(t)[,  ct(t)  =  h[et(t),  o(t)] 


(7) 


After  appropriate  differentiation  and  manipulation,  one  then  obtains 
the  desired  generalizations  of  eqs.  (6)  in  integral  form  as 


Thus,  the  axial  creep  strain  for  use  in  eq.  (4)  is  given  by  eq.  (5) 
with  eqs.  (8) . 

In  order  to  obtain  the  solutions  in  the  next  two  sections,  we 
shall  also  require  expressions  for  the  total  axial  strain  and  the 
total  lateral  strain.  For  the  total  axial  strain  e(t)  we  simply 
add  a  linear  elastic  term  to  the  axial  creep  strain,  i.e. 


c(t)  -  ~  +  E(.(t) 


(9) 


where  E  is  the  elastic  modulus.  And  for  the  total  lateral  strain 
f?(t)  we  follow  the  apporach  of  Courtine,  Cozzarelli  and  Shaw  [11], 
and  describe  lateral  contraction  by  means  of 


if'—*  -~ 
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e^(t)  =  -[v  ^  +  vg  eg(t)  +  v t  bt(t)]  (10) 

where  v  is  the  usual  linear  elastic  Poisson's  ratio,  while  v 

s 

and  v  are  steady  creep  and  transient  creep  Poisson's  ratios  respec¬ 
tively.  For  metals,  is  close  to  1/2  (incompressible)  whereas  v 

is  between  v  and  vg  .  Note  that  we  have  not  employed  above  the  simpli 
f ied  approach  of  combining  elastic  and  transient  creep  strain  into  a 
single  nonlinear  instantaneous  term  (e.g.,  see  [4]). 


2 . 3  Rupture  Time 

The  presence  of  the  creep  strain  in  the  damage  relation  (4)  was 
motivated  by  the  experimental  observations  in  [5,7,8],  and  results  in 
a  coupling  of  the  damage  and  creep  strain  relations  which  is  necessary 
for  a  physically  reasonable  estimation  of  rupture  time.  We  have  not 
followed  Kachanov's  approach  [6]  of  including  the  damage  in  the  creep 
strain  relation,  since  it  was  felt  that  this  additional  coupling  would 
complicate  the  formulation  unnecessarily.  Therefore  the  present 
formulation  does  not  account  for  tertiary  creep  at  constant  stress, 
but  this  is  not  necessarily  a  restriction  if  for  example  one  defines 
the  rupture  time  as  that  time  at  which  rapid  tertiary  creep  and  rupture 
commence.  Furthermore,  the  significance  of  tertiary  creep  up  to  the 
point  of  actual  rupture  can  be  minimized  if  care  is  taken  to  base  the 
creep  model  on  constant  stress  data  (as  done  here)  rather  than  on 
constant  load  data. 

Specifically,  let  t  be  the  rupture  time  corresponding  to  a 
critical  value  of  damage  for  rupture  .  For  the  case  of  constant 

stress  eqs.  (1),  (5)  and  (6)  then  yield 


where  gives  the  rupture  time  in  implicit  form.  If  we  ignore  transient 
creep  (p  -►  °°)  we  may  obtain  t  explicitly  as 


.  not .  ( 

=,  •  ch-K'^] 


which  agrees  with  the  result  given  in  [5],  except  that  there  oq  was 
the  initial  stress  for  constant  load.  It  was  also  pointed  out  in  [5] 
that  the  positive  quantity  (not+y) /(a+6)  is  analagous  to  the  exponent 
of  stress  (generally  designated  as  v)  in  the  Kachanov  damage  law. 

In  the  next  two  sections  we  examine  several  more  complicated  tests. 
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3.  DAMAGE  SOLUTIONS  FOR  VARIOUS  PRESCRIBED  STRESSES  AND  STRAINS 


Some  insight  into  the  physical  nature  of  the  proposed  creep  model 
may  be  gained  by  obtaining  explicit  expressions  for  the  damage  D(t) 
with  various  prescribed  stress  and  strain  histories.  Thus,  in  this 
section  we  obtain  solutions  to  eq.  (4)  for  the  constant  stress  rate, 
step  stress  and  relaxation  tests,  and  also  make  some  observations  on 
the  associated  rupture  times. 


3.1  Constant  Stress  Rate  Test 

As  a  first  example,  consider  the  case  of  a  stress  increasing 
linearly  with  time  from  a  zero  value  at  t  =  0  ,  i.e. 

o(t)  =  k  t  H(t) ,  k  >  0  (13) 


where  k  is  the  constant  stress  rate  and  H(t)  is  the  unit  step  func¬ 
tion.  Restricting  q  to  odd  positive  integers,  eqs.  (5)  and  (8)  yield 
for  this  case  the  axial  creep  strain 


,  ,  r>"  t”+1 .  ,k,q 

ec(t>  -  l<i>  “Sr +  (v> 


1  X 


11=0 


l-1)  T  3!  tq_i+  Tqq'e~t/T' 

(q-i) !  q  j 


H(t) 

(14) 


Eq.  (4)  with  eq.  (13)  then  yields  the  damage  in  terms  of  e^Ct) 


as 


■-P-  (U_ 

Cec(t)a 


kY  tY+6 


H(t) 


(15) 


If  we  ignore  transient  creep,  eqs.  (14)  and  (15)  yield  the  simple 
expression  for  damage 


i»+*rw'  w . 
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D(t)  =  -L-_  (-jL)  A-na  kY+nu  J+«(n+1>+5  H(t) 
(n+l)a  Y+6 


Note  that  in  eqs.  (15)  and  (16)  D(t)  -*■  ®  as  t  ->  «  ,  and  as  in  the 
constant  stress  test  there  will  always  exist  a  finite  rupture  time  t ^ 
at  which  the  damage  reaches  a  critical  value  . 


3 . 2  Step  Stress  Test 

Now  consider  the  case  of  a  constant  stress  applied  at  t  =  0  and 
then  suddenly  changed  to  K  0Q  at  t  =  t^  ,  i.e. 


o(t)  =  oo  {  H 


(t)  -  H(t  -t;)  +  KH(t  -  tx)y  ,  K  >  0  (17) 


In  this  case,  eqs.  (5)  and  (8)  give 


a  n 

e  (t)  =  ( r~)  [ tH(t)  +  (Kn  -  1) (t  -  t.)H(t  -  t.)] 

C  A  11 


+  (— ^) q  [  ( 1  -  e  t/l)H(t)  +  (Kq-l)(l  -  e  (t  tl)/T)H(t  -  t^ 


Next,  eq.  (14)  gives 


— ---  =  0QY[t*H(t)  +  (Kr-l)(t6-  t16)H(t  -  t]) 


Cec(t) 


Again,  if  we  ignore  transient  creep  a  simplified  expression  for  damage 
is  obtained  as 


^  k  '  ”  . i  „  ,  0--  *"  ■  ■  -  i -prf - >■  r*,' 
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D(t)  =  C  x"an  ooY+najt6+a  H(t)  +  [ (KY-1)  ta(t6-  0 

+(Kan-l)  t6(t“-  t1U)+(KY-l)(Kan-l)(tfi-t16)(ta-t1a)]H(t-t1)| 

(20) 

As  in  the  constant  stress  rate  test,  D(t)  -*■  °°  as  t  ->  °°  in  the 
present  case  and  a  finite  rupture  time  t  again  exists.  Also  note 
that  although  in  accordance  with  eq .  (17)  o(t)  is  discontinuous  at 
t  =  t^  ,  the  damage  D(t)  in  accordance  with  eqs.  (19)  and  (18)  is 
continuous  at  this  point  while  D(t^)  is  discontinuous . 

3 . 3  Relaxation  Test 

As  a  third  example,  consider  the  case  of  a  constant  strain 
suddenly  applied  at  t  =  0  ,  i.e. 


e(t)  =  co  H(t) 


(21) 


When  e(t)  is  prescribed,  eq.  (9) (with  eqs.  (5)  and  (8))  is  a  nonlinear 
differential  equation  in  o(t)  which  may  in  general  be  difficult  to 
solve.  However,  if  we  ignore  transient  creep,  this  differential 
equations  with  condition  (21)  simplifies  to 


I  £♦<!>  -» 


t  >  0 


(22) 


where  o(0+)  =  E  .  For  n  >  1,  the  solution  to  eq.  (22)  is  readily 
obtained  as 


o(t)  =  E 


-n+1 


+  (n-1) 


EnA_n  t 


-1  -n+1 


H(t) 


(23) 
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The  relaxation  test  with  transient  creep  neglected  is  thus  equiv¬ 
alent  to  a  test  where  stress  is  prescribed  in  accordance  with  eq.  (23) 
Substituting  eqs.  (21)  and  (23)  into  eq.  (9)  we  then  obtain  the  creep 
strain  as 


ec(t)  =  |co  -  [eo  n+1+  (n-l)En  A  n  t ]  i  H(t) 


Finally,  eq.  (4)  with  eq .  (23)  yields  the  damage  in  terms  of  tc(t) 


Cc  (t)' 
c 


-1  n  Vsf  6(n-l)-Y 

E  An  6r,,„  -  -n+1 . ,  , ,  -n. , 6(n-l) 

6(n-l)-Y J  |[(Ee0)  +(n-l)EA  t] 

6(n-l) — 6 

[(Eto)'n+l]6(n_1)  H(t) 


In  the  limit  as  t  -*■  «°  we  obtain  from  eqs.  (24)  and  (25) 


.  •  ^  a+Y-6(n-l)„Y-|5n/  An<5  1 

Lim  D(t)  Ce  E'  — 77 — rr  ,  w  >  0 

=  .  0  yr~& (n-1) J 

t-H»  CO  ,  (o  <  0 


with  u)  =  y  -6(n-l) 


We  thus  obtain  the  interesting  result  that  whereas  for  the  material 
power  combination  w  >  0  the  damage  may  not  necessarily  ever  obtain 
a  critical  value  ,  the  existence  of  a  finite  rupture  time  t  is 

guaranteed  when  to  <  0.  A  similar  result  is  obtained  when  one  studies 
damage  during  stress  relaxation  using  the  Kachanov  damage  formulation. 

We  treat  the  constant  load  test  separately  in  the  next  section, 
since  that  case  is  somewhat  more  complicated  than  the  three  cases 
considered  in  this  section. 
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4.  DAMAGE  SOLUTIONS  FOR  CONSTANT  LOAD 


Consider  now  a  bar  of  instantaneous  cross-sectional  area  A(t) 
under  a  constant  load  PQ  suddenly  applied  at  t  =  0.  The  bar 
experiences  a  continuously  increasing  stress 


o(t) 


A(t) 


(27) 


due  to  the  effect  of  lateral  contraction,  and  the  existance  of  a  finite 
rupture  time  is  assured.  Employing  the  logarithmic  definition  for  large 
strain,  we  write  the  lateral  strain  rate  as 


•  ,„x  _  1  A(t) 

"  R(t)  2  A(t) 


(28) 


where  R(t)  is  the  instantaneous  cross-sectional  radius.  In  this  sec¬ 
tion  we  determine  the  solutions  for  A(t)  ,  o(t),  &c(t)  anc^  D(t) , 

first  for  the  case  of  large  deformation  with  transient  creep  excluded 
(or  included  approximately  in  the  elastic  term)  and  then  for  the  case  of 
small  deformation  with  transient  creep  included. 


4.1  Large  Deformation  -  Transient  Creep  Excluded  (or  Included 
Approximately) 

For  this  case  we  set  =  0  in  eq.  (10)  and  also  assume  that 
the  steady  creep  is  incompressible  (vg  =  1/2)  ,  whereupon  with  the  use  of 
eqs .  (8)  and  (27)  we  obtain 


:,(t)  =  -v 


1 

2 


Po  X 
A(t) A  j 


t  >  0 


(29) 


This  result  in  combination  with  eq.  (28)  leads  to  the  differential  equa¬ 
tion  in  A(t) 


2vP 


°  .n-2 
—  A  dA 


An-1dA,  t  >  0 


(30) 


One  may  obtain  a  rough  correction  for  the  exclusion  of  transient  creep 

•  ic 

by  replacing  E  with  E  ,  a  smaller  fictitious  modulus  of  elasticity. 
The  initial  condition  for  eq.  (30)  is  prescribed  at  t  =  0+  (i.e., 
immediately  after  the  load  is  applied) ,  i.e. 


A(0+)  =  A 

o 


(31) 


Since  the  strain  at  t  =  0+  is  small  we  may  write 
tj,  (0+)  -  1/2[A(0+)  -  A(0“)]/A(0+)  and  thereby  obtain 

2vP 

A(0+)  *  A(0  )  -  (32) 


from  which  we  see  that  A(0+)  differs  from  the  known  A(0  )  by  a 
very  small  amount . 

Differential  equation  (30)  with  initial  condition  (31)  may  be 
integrated  to  yield 


2nvo 


o 

(n-l)E 


1 


2nvo 


(n-l)E 


o  n 

n(y^)  t,  t  >  0 


(33) 


where  =  Pq/Ao  .  Eq.  (33)  gives  the  function  AQ/A(t)  =  F(t)  in 
implicit  form,  but  explicit  expressions  may  be  obtained  for  approxi¬ 
mate  cases.  For  example,  if  we  neglect  the  elastic  term  (E  -*•  °°) 


•> — — n, - - 
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eq.  (33)  gives  the  simple  relation 


A 


o 

A(t) 


'V 


1  -  n(— )  t 


1 

n 


t  >  0 


(34) 


On  the  other  hand,  we  can  retain  the  elastic  term,  but  approximate  Aq/A 
within  the  brackets  by  unity,  and  accordingly  obtain 


1 


o  — 

n(n-l)E  ,_o, 

(n-l)E-2nvo  a 
o 


1 

n 


t  >  0 


(35) 


Since  eq.  (35)  contains  eq.  (34)  as  a  special  case,  we  shall  use  eq. 
(35)  in  the  remainder  of  this  subsection,  with  E  replaced  by  E* 
when  transient  creep  is  approximated. 

The  stress  now  follows  immediately  from  eq.  (35)  as 


o(t)  =  o 


o  A(t) 


&  a 


1  - 


0  n  - 

n 


n(n-l)E  ,  o, 
(n-l)E-2nvo  '‘A  ’ 


t  >  0  (36) 


Next,  the  creep  strain  as  obtained  by  substituting  eq.  (36)  into  the 
first  of  eqs.  (8)  is  given  by 


e 

c 


(t) 


-In 


1- 


n(n-l)E 


(n-l)E-2nvo  'A 
o 


o  n  -i 

(t^)  t 


n(n-l)E 


(n-l)E-2nvo 


t  >  0 


(37) 


Note  that  e  (t)  -►  <*>  as  t-+t  ,  where 
c  max 


t 

max 


(n-l)E-2nvo  o  _1 

- °  (_°) 

n(n-l)E  A 


(38) 
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which  is  analagous  to  the  ductile  rupture  criterion  due  to  Hoff  112]. 

Finally,  the  damage  is  obtained  in  terms  of  c^(t)  from  eqs.  (4) 
and  (36)  as 


„°iii  %  r  < 

Ce  (t)  ‘ 

c 


,  6n  , 

( — r~)  t 

Y—on  niaA  jji 


ax  i  -  r~) 
u  max 


tof 

6n 


t  >  0 


(39) 


J 


Note  that  as  expected  a  finite  rupture  time  tf  exists  corresponding 

to  a  given  critical  damage  ,  but  now  in  contrast  with  the  constant 

stress  test  [eqs.  (11-12)]  D(t)  -*■  «>  as  t  -*■  t  .  Also  note  that  in 

max 

general  y  -  6 n  in  eq .  (39),  and  in  fact  for  the  very  special  case 
y  =  6n  eq .  (39)  simplifies  to 


D(t) 

Cec(t) 


a  '  o 


ln(  1  -  -r—) 1 


t 

max  J 


(40) 


In  order  to  illustrate  the  character  of  the  above  strain  and  damage 

solutions,  we  employ  the  data  yiven  in  [8]  for  the  following  two  steels: 

material  a:  AISI  310  stainless  steel  at  600°C 

material  b:  2.25Cr  IMo  ferritic  steel  at  550°C 

The  various  material  parameters  (elastic  -  E*,v;  creep  -  n,X  n;  damage 

-a,Y,<$,C)  for  these  two  materials  are  summarized  in  Table  1,  where  the 

E*  indicates  we  are  using  the  rough  correction  for  the  exclusion  of 

transient  creep.  Using  these  values  with  an  initial  stress  cq  equal 
2  o 

to  30  kg/mm  for  material  a  and  25  kg/mmz  for  material  b,  we  obtained 

the  constant  load  curves  r.  (t)p  [eq.  (37)]  and  D(t)p  [eq.  (39)] 

c  o  o 

shown  plotted  in  Figures  1  and  2  respectively  vs.  time  in  hrs.  For 

purposes  of  comparison  we  have  also  shown  in  Figures  1  and  2  the  constant 

stress  curves  e  (t)  and  D(t)  [eqs.  (1,5,6)  with  e  =  0]  for  these 
c  o  o  ’  t 
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same  values  of  the  material  constants  and  o  .  It  is  noteworthy  that 

o  1 

the  curves  at  constant  load  differ  very  significantly  from  those  at 
constant  stress. 

In  Figure  L  we  have  indicated  the  ductile  rupture  times  (t  ) 

°  max 

corresponding  to  <  •  ,  Also,  in  Figure  2  we  have  marked  the  rupture 

times  at  constant  load  and  stress  [^tr)p  >  (tr)0  1  corresponding  to 

°  -4° 

critical  values  of  damage  Dr  equal  to  25  x  10  for  material  a  and 

_4 

15  x  10  for  material  b  (from  (8]).  These  various  rupture  times  are 
summarized  in  Table  2;  as  expected  the  rupture  times  at  constant  load 
are  less  than  both  the  ductile  rupture  times  and  the  rupture  times  at 
constant  stress. 


TABLE  1  -  Material  Constants 


/ 


17 


i 


TABLE  2  -  Rupture  Times 

Constant  Load  and  Stress  -  Large  Deformation  - 
Transient  Creep  Included  Approximately  in  E* 


[ 

Material 

a 

b 

J 

t  ,  hr 

max 

86.4 

57.8 

ur)P  .  hr 

o 

77.6 

33.6 

(tr>„  '  hr 
o 

196.1 

55.9 

4 . 2  Small  Deformation  -  Including  Transient  Creep  as  a  Separate  Term 

Consider  again  the  derivation  of  a  differential  equation  in  A(t)  , 

except  that  now  we  retain  the  transient  creep  term  et(t)  and  also  do 

not  restrict  v  to  1/2.  To  this  end  we  use  eq.  (28)  to  form  the 
s 

combination 


••  ,  .  ,  1  *  .  .  1  ,A.  1  A 

+  7  eS-(t)  =  2 


2t  A 


(41) 


If  we  now  substitute  eq.  (10)  with  eqs.  (8)  and  (27)  into  the  left  hand 
side  of  eq.  (41),  we  get  the  nonlinear  differential  equation 


vP  vP  P  n  —  v  P  n  v_  P  9  — 

-r  8  +  *r  s  +  <r>  “n+r<r>  “°  +  r<r)  B 


■  i  <!>  +  it  <!> 


t  >  0 


(42) 
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where  B(t)  =  l/Mt)  • 

Equation  (42)  cannot  in  general  be  solved  in  closed  form,  and  thus 
we  shall  simplify  the  problem  by  introducing  the  assumption  of  small 
deformation.  Accordingly,  we  write 


A(t)  =  Aq  -  A(t)  ,  A(t)  «  Aq 


(43) 


where  A(t)  is  a  small  area  increment  and  Aq  has  been  defined  by 
eq.  (31)  with  eq.  (32).  Substituting  eq.  (43)  into  eq.  (42)  and  neg¬ 
lecting  higher  order  terms,  we  now  obtain  a  linear  differential  equa¬ 
tion  in  A(t)  as 


,  vo  .. 

<i "  ir)A  + 


2t 


VOo  f°o"  qvt  ,°o. 

TT  -  ""s  <->  -  —  (r> 


nv 


a  n  A  v  a  n 

,  O.  r  O  S  ,  O. 

(__)  A  _  _  (_)  , 


t  >  0 


(44) 


where  as  before  oo  =  Pq/Ao  • 

For  convenience  we  introduce  the  nondimens ional  variables 


o 


Equation  (44)  then  becomes 

a^  A"  +  a2  A'  -  a^  A  =  a^/n  ,  t  >  0  (46) 

where  the  prime  indicates  d/dt  ,  and  the  nondimensional  coefficients 
are  defined  by 


— *  , 


*3^.. 
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i 


1 

2 


vo 

o 

E 


(47a) 


on  o  9 

T(^)  -  qvt0 


(47b) 


a  n 
,  °\ 

a  =  nv  t(-t— ) 

3  s  A 


(47c) 


As  initial  conditions  for  eq.  (46)  we  require  A(0+)  and  A'(0+)  . 
The  former  tollows  directly  from  eqs.  (31),  (43)  and  (45)  as 


A(0+)  =  0 


(48) 


In  obtaining  the  latter  we  first  note  that  eg(0+)  -  et(0+)  =  0  ,  and 
use  eqs#  (8),  (10)  and  (28)  to  obtain 


4<0+) 


1  Mp+)  /  o. 

2  A—  ~  W 

o  o 


A(0+) 


o  n 
(  °\ 
Vr0 


(49) 


Introducing  eqs.  (43),  (45)  and  (47)  into  this  result  and  again  assuming 
small  deformation  we  obtain  the  desired  condition  as 


*'<0+>  =  !~n  +  \ 


V 


a2_ 

alq 


(50) 


The  solution  to  differential  equation  (46)  with  initial  conditions 
(48)  and  (50)  is  easily  obtained.  For  example,  in  the  special  case  of 
q  =  n  we  get 


A  ( t ) 


1 

n 


I  -a„+  i+  [ 
4  \  2nA 


+ 


/v  2V  A 

l  ~  2nA 


t  >  0 


(51) 
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where 


J 


v. 


A  =  (ai:  +  4 


ala3} 


(52) 


The  stress  then  follows  directly  from  eq.  (51)  in  accordance  with 


o(t)  =  oQ[l  +  A(t) ]  ,  t  >  0 


(53) 


Substituting  eq.  (53)  into  eq.  (5)  with  eq.  (8)  and  neglecting 
higher  order  terms,  we  obtain  the  creep  strain  as 


ec(t) 


a,  on 
1  .  o. 

r3~  '<-> 


-a„+  2a, +  A  r  2a 


““i' 

+  (-^2Z— > 


n  2a, 


t(-A-t)(t£)  +  c 


a  °-!  -a„  +  A 


-a2+  A  'X 


)(~) 


2a  -  a 2+  A  y 


exP(— —  t) 


a„-  2a, +  A  r  2a,  o  n  2a 


2  1 
+  <■  > 


o 


tC-1— -)(-"-)  +  (y- - - - r)(— ) 

~a2~  ^  ^  2a^-  a^~  A  p 


“a2"  A 

exp(~2ir  E) 


t  >  0 


(54) 


Although  the  creep  model  at  constant  stress  (i.e.,  eqs.  (5)  and  (6))  does 
not  contain  a  tertiary  creep  stage,  eq.  (54)  does  yield  a  region  of  increas¬ 
ing  creep  rate  after  a  critical  time  t  .  Setting  ec"(tc)  =  0  we 
obtain  this  critical  time  as 


t 

c 


1  n< 


0  n  on 

(-a2~  A)[(2aj-  a2~  A)x(y^)  +(-a2~  A)  (-^) 


n 


(~a2+  A)  [  (2a^-  a2+  A)x(— )  +(-a2+  A)  . 


n  [ 

0 

(— )  1  ’ 


(55) 
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Finally,  the  damage  for  small  deformations  is  obtained  from  eq.  (4) 
as 


1 


„  /-v  01 

Ce  (t) 

C 


-£  -i 

6n  6n 


-a0+  2an  +  A  2a.  -a0+  A 

+  2  ,  ?  -)(-v-r)«.p<-|— >t 

6  2nA  -a2+  A  2a^ 


a0-  2a,+  A  2a 

+  H(-2 - i - )  ( 

<5  ^  2nA  M 


~a_-  A 


1  ,  ,  2  ,-1 

- -)  exp  ( — x~— )  t  > 

~a2~  A  2al  J 


(56) 


where  we  again  find  that  t  is  finite. 

In  order  to  illustrate  the  above  solutions  the  creep  strain  was 

calculated  using  the  parameters  assembled  in  [13]  from  the  data  of  Garafalo 

et  al.  [9,14]  for  the  following  steel: 

material  d:  AISI  316  stainless  steel  at  816°C 

These  elastic,  steady  creep  and  transient  creep  parameters  (E,n,A  n,q, 

p  q,r)  are  given  in  Table  3  along  with  the  additional  assumed  Poisson's 

ratios  v  =  .3,  =  .4  and  vg  =  .5  .  Figure  3  shows  the  e^  curves  for 

constant  load  [eq.  (54)]  and  constant  stress  [eqs.  (5,6)]  for  equal 

2 

to  7.03  kg/mm  ,  and  we  again  see  that  these  two  solutions  differ  by  a 

considerable  amount.  Upon  comparing  Figures  1  and  3  we  see  that  the 

ductile  rupture  time  t  is  no  longer  finite  when  deformations  are 

max 

small.  However,  we  now  have  the  critical  time  t£  [eq.  (55)]  corres¬ 
ponding  to  the  point  of  inflection,  and  in  the  present  example  this 

occurred  at  84  hours  (see  Figure  3) . 


X' 
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TABLE  3  -  Material  Constants 


Material  d 

2 

E,  kg/mm 

10,500 

.3 

V 

.4 

t 

V 

.5 

s 

n  =  q 

3.5 

A  n,  mm^n/kgn-hr 

0.795  x  10~6 

p  q,  mm2q/kgq 

0.802  x  10~4 

T  ,  hr 

21.7 

5.  DISCUSSION 

Starting  with  a  damage  relation  for  constant  one-dimensional  tensile 
load  given  in  [5],  an  extended  damage  relation  was  proposed  for  time 
dependent  one-dimensional  tensile  stress  histories.  This  relation  con¬ 
tains  a  power  of  the  strain  history  and  also  a  power  of  a  time  integral 
over  a  power  of  the  stress  history,  thereby  ensuring  a  continuous  damage 
history  even  for  a  discontinuous  stress  history.  Employing  an  axial 
strain  relation  containing  both  steady  and  transient  creep  terms  with  a 
Lateral  strain  relation  containing  three  Poisson's  ratios,  damage  solu¬ 
tions  were  obtained  first  for  the  constant  stress  rate,  step  stress  and 
relaxation  tests  and  then  for  the  constant  load  test. 

In  all  tests  considered,  except  the  relaxation  test,  there  existed 
a  finite  rupture  time  corresponding  to  a  critical  value  of  damage;  the 
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existence  oi  a  finite  rupture  time  was  guaranteed  in  the  case  of  the 
relaxation  test  only  when  a  material  power  combination  was  negative. 
Constant  load  test  solutions  were  obtained  both  for  the  case  of  large 
deformation  with  transient  creep  excluded  (or  included  approximately 
in  the  elastic  terms)  and  for  the  case  of  small  deformation  with  tran¬ 
sient  creep  included  as  a  separate  term;  in  both  cases  the  creep  strain 
and  damage  solutions  are  considerably  greater  than  the  corresponding 
solutions  at  constant  stress.  in  the  former  case  a  ductile  rupture  time 
was  obtained  corresponding  to  infinite  strain.  In  the  latter  case  a 
critical  time  was  observed  marking  the  beginning  of  a  region  of  increas¬ 
ing  creep  rate. 

A  testing  program  for  the  determination  of  the  material  parameters 
would  most  logically  begin  with  the  one-dimensional  constant  tensile 
stress  test,  with  the  axial  creep  strain  and  damage  measured  as  a  func¬ 
tion  of  time.  Density  variation  is  one  useful  index  of  damage,  but  it 
may  be  necessary  to  also  use  other  methods  such  as  the  actual  counting 
of  voids.  The  most  reliable  procedure  would  probably  be  to  first  fit 
the  axial  creep  strain  data  to  determine  the  material  parameters  in  the 
model  for  and  then  to  smooth  the  damage  data  with  the  calculated 

> c(t)  before  finally  fitting  this  damage  data  to  determine  the  material 
parameters  in  the  model  for  D(t)  .  Then  the  validity  of  the  proposed 
model  for  one-dimensional  time  dependent  tensile  stress  could  be  checked 
by  performing  one  or  more  of  the  tests  discussed  in  this  paper,  with 
the  step  stress  and  constant  load  tests  being  possibly  the  most  useful. 
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APPENDIX:  LIST  OF  SYMBOLS 


,  cl 

A(t) 

A 

o 


A(t) 

A(t) 

B(t) 

C 

D(t) 

D 

r 

E 


E* 

H(t) 

k 

K 


n,q 


P 

o 

R(t) 

t 

t 

max 


t 

r 


t 


t 

c 


nondimensional  coefficients,  eqs.  (47) 
instantaneous  circular  cross-sectional  area 

initial  cross-sectional  area  immediately  after  load  application, 
eq.  (31) 

small  area  increment,  eq.  (43) 
nondimensional  area,  eq.  (45) 
reciprocal  of  A(t),  eq.  (42) 
damage  constant,  eq.  (1) 
damage,  eq.  (4) 

critical  value  of  damage  for  rupture 
elastic  modulus 

fictitious  elastic  modulus  to  approximate  transient  creep 

unit  step  function 

constant  stress  rate,  eq.  (13) 

constant  multiplier  in  step  stress  test,  eq.  (17) 

creep  constants,  eq.  (6) 

constant  load  suddenly  applied  at  t  =  0 

instantaneous  cross-sectional  radius 

time 

ductile  rupture  time,  eq.  (38) 
rupture  time 

constant  time  in  step  stress  test,  eq.  (17) 
nondimensional  time,  eq .  (45) 

critical  time  for  constant  load  at  beginning  of  period  of 
increasing  creep  rate,  eq.  (55) 
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i 


e(t) 


ec(t) 

tg(t) ,£t(t) 


h(t) 


O 

A,p,t 

v 

VVt 

p(t) 

Po 
a(  t) 


o 

0) 

(•) 

(  )' 

(  )p 

o 

(  >o 

( i°( ) 
(  ) 


b’ 

d 


damage  constants,  eq.  (1) 
nondimens ional  constant,  eq.  (52) 
total  axial  strain,  eq.  (9) 
axial  creep  strain,  eq.  (5) 

axial  strains  due  to  steady  creep  and  transient  creep 

respectively,  eqs.  (6) 

total  lateral  strain,  eq.  (10) 

constant  strain  in  relaxation  test,  eq.  (21) 

creep  constants,  eqs.  (6) 

elastic  Poisson's  ratio 

steadt  creep  and  transient  creep  Poisson's  ratios,  eqs.  (10) 
density 

initial  density 

one-dimensional  tensile  stress 
constant  stress 

material  power  combination  in  relaxation  test,  eq.  (26b) 
indicates  d/dt 
indicates  d/dt 

indicates  constant  load  test 
indicates  constant  stress  test 

indicates  materials  a,b  and  d 
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Key  to  Subscri p_ts : 

a  -AISI  310  (600°C) 
b  -2.25  Cr  IMo  (55( 


Fig  1  -  Strains  at  constant  load  and  stress  -  Large  deformation 
-Transient  creep  included  approximately  in  elastic  term 


Key  to  Subscripts : 

a  -AISI  310  (600°C) 
b  - 2.25Cr  IMo  (550°C) 


Fig  2  -  Damage  at  constant  load  and  stress  -  Large  deformation 
-  Transient  creep  included  approximately  in  elastic  term 


Key  to  Subscripts 
d-AISI  316  (816°C) 


Fig.  3-  Strains  at  constant  load  and  stress -Small  deformation 
-Transient  creep  included  as  separate  term 
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